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Considering a double-barrier structure formed by a silicon quantum dot covered by natural oxide with two 
metallic terminals, we derive simple conditions for a step-like voltage-current curve. Due to standard chemical 
properties, doping phosphorus atoms located in a certain domain of the dot form geometrically parallel current 
channels. The height of the current step typically equals to N x 1.2 pA, where N = 0, 1, 2, 3 . . . is the number 
of doping atoms inside the domain, and only negligibly depends on the actual position of the dopants. The 
found conditions are feasible in experimentally available structures. 
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The fabrication of Si nanostructures became possi- 
ble through recently developed new technologies [Q, ||. 
Individual silicon quantum dots (SQD) reported in 
H are spherical Si particles with diameters d in the 
range 5-12 nm covered by a 1 - 2 nm-thick natural S1O2 
film. Metallic current terminals made from degenerately 
doped Si are defined lithographically to touch each in- 
dividual dot from above and from below. 

To ensure metallic electrodes the donor concentra- 
tion n should be n > riMott, where riMott = 7.3 x 
|10 17 cm~ 3 . The critical concentration nMott is defined 
■by the Mott criterion [Q], introducing the transition to 
a metallic type of conductivity in a semiconductor at: 



a B x (n M ott) 1/3 = 0.27. 



(1) 



where a B is the Bohr radius of an electron bound to a 
donor inside the Si crystal, in the case of phosphorus- 
donors a B = 3 nm H . 

As for the doping of the dot, the situation concerning 
a Mott transition in that small dots is much less trivial 
than the one described by Eq. (|]) . Let us consider dots 
with diameters d = 10 nm formed from n-doped Si with 
n = riMott as an illustrative example. Then each dot 
contains in average one donor. Note that we will con- 
sider degenerately n + -doped electrodes with n 3> riMott 
which ensures metallic conduction up to the borders of 
the dot. 

Real fabrication technology provides a wafer with 
hundreds of SQDs on it with current leads towards each 
individual SQD. Dots in average have the same value of 
mean dopant concentration n, which is determined by 
the parent material of bulk silicon the dots are formed 
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from. However, on the level of each individual SQD 
we will always have exactly integer number of doping 
atoms. If, as in the example above, the average num- 
ber of dopants iV tot = 1 the actual number of donors in 
the dot can have values N to t = 0,1, 2, 3,..., with values 
larger than these very unlikely. 

Our objective is to illustrate, that SQDs from the 
same wafer fall into several distinct sets of approxi- 
mately the same conductance. The typical value of con- 
ductance for each set is nearly completely determined by 
the number TV of donors present in a certain part of a 
SQD so that N labels each set of SQDs. 

Summarizing the above, we need for a quantization 
of the conduction through a dot with N donors the fol- 
lowing conditions: 

• Size d of the dot comparable with Bohr radius: 
2 < d/a B < 5. 

• Average doping n of the dot n < d~ 3 , leading 
to a mean number of dopants Ntot < 1, so that 
Ntot = 0, 1, 2 are the most probable configurations 
of an individual SQD. 

• Doping of the electrodes n e i ^> riMott, so that cur- 
rent leads are perfectly metallic. 

• Dot covered by an oxide layer thick enough to sup- 
press ballistic transport through the dot. 

In fact all these condition can be simultaneously sat- 
isfied for SQD fabricated with the method mentioned 
above ||. 

Model system. We use a simple model of a cu- 
bic SQD with d > 2a B (we will use d = 10 nm for 
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estimates), covered with an oxide layer with thickness 
S = 2nm, height || B = 3.15 eV and contacted with 
current terminals from left and from right. The x-axis 
is oriented from left to right along the current flow, as 
shown in Fig.|l|. 
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Figjl]. Potential profile of the dot covered by an oxide 
layer at V = 0. A donor is marked with a short bar 



A tunneling current is injected into the dot via the 
oxide barrier from the top (source at x = 0) and leaves 
the dot at the bottom (drain at x = d). Due to the pres- 
ence of the oxide barriers this current is non-ballistic 
and non-thermal. We assume that the high potential 
barriers associated with the oxide layers are not much 
affected by the voltage and the tunneling charges. We 
concentrate on what happens between these effective 
source and drain (Fig.^), as in ||. 

In the case when the dot can be regarded as an insu- 
lating system it is reasonable to assume that the applied 
voltage equally drops over the potential barriers and the 
dots. For simplicity we neglect the difference of the di- 
electric constants of the oxide barriers and the dot. In 
this approximation we can introduce an effective volt- 
age V e ff = V(d — 28) I d — 0.6 V describing the part of 
the total transport voltage V applied between effective 
source and drain which drops across the dot itself. 

In this rude approximation we neglect the effect of 
spatial quantization upon values on the ionization en- 
ergy, the conductivity gap and material parameters of 
silicon. 

Dot without donors. At V e ff = the Fermi level 
inside the dot is situated in the middle of the gap, i.e. 




Fig.g. Potential profile of the dot between effective 
source and drain biased with Veff = V\ (thick solid 
line) and V e f / = Vth (dashed line) 



Eg/ 2 bellow the conduction band edge (E g = 1.14 eV 
at 300 K). 

As V e ff grows, the bottom of the (still empty) con- 
duction band bends down accordingly. When the con- 
duction band in the dot close to the drain aligns with the 
Fermi level of the emitter we expect a drastic increase 
in the tunneling current. This threshold Vth voltage 
(Fig.|) for V ef f is given by V th = E g /(2e), regardless of 
the number N tot of dopants in the dot (as long as the 
dot is not yet metallic, of course). In the following we 
therefore limit our studies to voltages 
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< V th = E g /(2e) = 0.57V. 



(2) 



In this voltage range we have a d-thick barrier 
(formed by the dot) with always finite height between 
effective source and drain. The intrinsic concentration 
of electrons and holes at 300 K is 1.4 x 10 10 cm -3 . Even 
at this high temperature the probability to have at least 
one intrinsic electron in a dot with size d = 10 nm is only 
1.4 x 10 -8 . So we would expect virtually no current in 
this mode. This is confirmed by direct electrical tests 
of SQD with the required properties. 

Single-donor channel. Let us now consider one 
single donor in the dot located at x with ionization en- 
ergy H Ed — 0.045 eV (for P as a donor). 

At zero temperature current is due to resonant tun- 
neling via non-ionized donor (as in Q for example). Dif- 
ferential conductance g(e) for the states with energy s 
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where r^ r is linewidth of Is state of electron bound to 
the donor due to coherent mixing with conduction states 
to the left (right) of the left (right) tunnel barrier. 

Oxide barriers (with height B — 3.15 eV and width 
6 = 2 nm) give dominant contribution to Ti_ r compared 
to contribution of the body of the dot (with typical 
height < E g /2 = 0.57eV and width < d = lOnm). So, 
we can approximate Ti_ r with linewidth T for the case 
of an impurity localized at distance 5 inside rectangular 
one-dimensional tunnel burrier kJ of height B: 



T r = r 



p% 



2pfk exp(—2n5) 

TT-B- 



kS 



2.5 x 10~ 9 eV, 
(4) 

where m is (true, not effective) electron mass, n = 

(2mB) 1 / 2 /fi, and p F = {2mE F ) 1/2 /U = (3n 2 n el ) 1/3 is 
Fermi wave number in the contact electrodes. The nu- 
merical estimate in (|^) is given for the electrodes doped 
up to n e i = 10 21 cm" 3 as in B|. 

Within approximation (0) point e = Sd brings func- 
tion g(e) given by (|^) to a sharp maximum g(sd) = 
e 2 /Trh of width T <C V t h. 

From Fig.|^ it is clear that resonant energy Ed = 
Eg/2 — xV e f f I (2ed) — Ed- This means that as soon as 
effective Fermi level eV e ff/2 reaches a certain threshold 
eVi/2, tunnel current J flowing through the structure 
acquires a step-like increase of 



Ji = - 

e 



eV/2 ( \r\ 9(e d )27TT 2eT 

g(e)de = = — = 1.2 pA, 

-eV/2 e n 



(5) 

If the impurity is located near the drain, i.e. d—as < 
x < d (as donor 1 in Fig.||), then threshold V\ for the 
effective voltage V e f / is given by 



V x = E g /(2e) - E d /e = 0.525 V. 



(6) 



In contrast, for an impurity located at distances 
Ax > 2dEd/E g from the drain (i.e. further away than 
the threshold case of donor 2 in Fig.||), no additional 
current channel via a single impurity can be opened at 
low enough voltages defined in (||) where virtually no 
background current is present. In the present case this 
value Ax = 0.8 nm, which returns us to the above crite- 
rion: only impurities located in the immediate vicinity 
(defined within the accuracy as) of the drain contribute 
to the single-impurity channel. 

This shows that in first approximation the conduc- 
tance of this channel does not depend on x. As shown 
above, a single-impurity channel already only selects im- 
purities located within a very narrow range of x close 
to the drain. 



Two-, three-, multi-donor channel. The above 
consideration shows, that due to the bend of the bot- 
tom of conduction band following the transport voltage, 
there is no chance to notice current flowing through a 
sequential chain of impurities (such as donors 1 and 3 in 
Figj^), connecting source and drain. The contribution 
of such a chain will be totally masked by the current 
flowing directly via the conduction band. The only way 
for multiple impurities to manifest themselves in quan- 
tized conductance is to form multiple geometrically par- 
allel singe-impurity channels situated close enough to 
the drain as considered above. 

Therefore, if N > 1 impurities fall into the thin layer 
near the drain to approximately the same x coordinate 
as that of donor 1 in Fig.|| (within the Bohr radius), 
we will see a switching-on of an A^-fold channel with 
current 



J N = NJi = N2eT/h = N x 1.2 pA 



(7) 



at the same threshold voltage V e ff = V\ = 0.525 V as 
for a single-donor channel (Fig.||). 
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Fig.g. Current-voltage characteristics of a model sys- 
tem (not to scale) 

All the above considerations are only valid as long as 
the dot itself can be regarded as an insulating system. 
As the number of donors in a SQD grows, the dot be- 
comes a metallic particle, and the conduction band edge 
in the dot aligns with the Fermi level of the electrodes. 
In a very simple estimate we define this transition to 
a metal when the total volume of N to t donors with an 
individual volume of 47r/3 x a\ exceeds the volume of 
the dot. This is an exaggerated version of the Mott cri- 
terion (Q) which holds not only in bulk, but in a small 
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structure, too. For the analyzed example from above 
this gives N tot — 8 as a limiting value. The practically 
interesting set 0, 1, 2, 3, . . . for both Ntot and N consid- 
ered above is still far bellow this limit. 

Quite a number of other mechanisms of electron 
transport might take place in this system. Surprisingly 
even taking into account such other mechanisms j?| does 
not change much the main idea of the present paper. 

In small dots with diameter d < 2cib = 6 nm the 
domain with iV active dopants extends to the whole 
dot, and thus N — N to t- In large dots with diameter 
d > 2ag — 6 nm the domain with active dopants is less 
than the dot itself and is localized near the drain. Hence 
the position of the domain, number N and value Jn all 
can be different for the current flowing in different di- 
rections. Really, when sign of applied voltage changes, 
the source and the drain change places. 

In a certain sense the discrete increase of dot's con- 
ductivity which follow the increase of the dopants num- 
ber could be regarded as a mesoscopic analog of the Mott 
transition between insulating and conducting states of 
the system. 

Useful discussions with I. Devyatov, M. Kupriyanov, 
and S. Oda are gratefully acknowledged. 
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